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Abstract 


The  mixed  potential  electric  field  integral  equation  (MPIE)  method  is  a  widely  used  numer¬ 
ical  method  in  wave  scj-ttering  analysis.  The  MPIE  and  all  other  classical  Maxwell  equation 
scattering  methods  involving  the  scalar  and  vector  potentials  are  based  on  the  Lorentz  gauge. 
In  this  report  we  analyze  a  series  of  waveguiding  and  scattering  problems  formulated  by  en¬ 
forcing  the  Coulomb  gauge.  Coulomb  gauge  Green’s  functions  are  derived  for  a  rectagular 
waveguide,  free  space,  a  cylinder  and  a  wedge.  With  the  MPIE  method  the  currents  on  a 
scatterer  are  shown  to  be  gauge  independent.  Comparisons  are  therefore  made  between  the 
computation  times  necessary  for  calculating  the  current  on  a  straight  thin  wire  and  arbitrary 
scatterers  in  the  presence  of  a  circular  conducting  cylinder  and  conducting  wedge.  Conclu¬ 
sions  are  drawn  concerning  the  utility  of  the  Coulomb  gauge  as  compared  to  the  Lorentz 
and  other  gauges. 
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Chapter  1 
Introduction 


Computers  have  allowed  researchers  to  solve  for  the  electromagnetic  fields  radiated  and 
scattered  from  structures  with  greater  physical  detail  than  could  be  accounted  for  by  purely 
analytical  means.  However,  the  variety  of  anaJytical  formulations  developed  under  various 
geometrical  constraints  as  well  as  computer  storage  and  run  time  limitations,  have  led  to 
a  proliferation  of  numerical  solution  methods.  No  single  numerical  method  has  been  found 
that  can  be  efficiently  applied  in  all  situations,  so  while  certain  techniques  are  widely  used, 
research  continues  on  strategies  to  improve  analytical/numerical  formulations  which  will 
shorten  computation  time  and  extend  the  physical  sophistication  of  the  scattering  body. 

One  numerical  method  that  has  gained  general  acceptance  is  the  integral  equation  ap¬ 
proach  which  models  conducting  or  dielectric  bodies  by  surface  patches  and  by  volumetric 
cubes  or  tetrahedrons.  Integral  equation  models  are  designed  to  compute  currents  from 
which  impedances  and  near  and  far  fields  can  be  calculated.  The  most  widely  used  integral 
equation  numerical  technique,  that  is  especially  useful  in  near  field  scattering  analysis  but 
compariable  to  other  methods  when  the  far  field  or  radar  cross  section  is  to  be  calculated, 
is  the  mixed  potential  electric  field  integral  equation  (MPIE)  method  described  for  example 
in  Rao  et.al.  [1].  The  MPIE  is  formulated  by  expressing  the  electric  and  magnetic  fields  in 
terms  of  the  electric  vector  potential  (A)  and  scalar  potential  (<^)  which  are  composed  of  in¬ 
tegral  and  differential  operators,  the  current,  and  the  free  space  Green’s  function.  Boundary 
condition  are  enforced  resulting  in  a  integral  equation  in  terms  of  an  unknown  volume  cur¬ 
rent  or  a  surface  current,  located  on  the  boundary  between  dielectric  regions  or  on  conductor 
surfaces,  that  is  solved  for  by  the  method  of  moments  [1]. 
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The  disadvantages  of  the  MPIE  method  are  that  the  scalar  potential  operator  contains 
two  derivatives  and  both  the  scalar  and  vector  potentials  have  a  singular  kernel.  These  are 
overcome  by  careful  selection  of  the  basis  functions  that  represent  the  current  and  charge  and 
the  testing  functions  that  produce  a  linear  set  of  solvable  equations  and  standard  numerical 
methods  for  singularity  removal.  The  advantages  of  the  MPIE  method  are  that  it  can  be 
applied  to  both  open  and  closed  bodies  and  the  basis  function  subdomain  dimensions  can 
be  adjusted  to  give  more  detail  in  regions  where  the  current  changes  rapidly  with  spacial 
position. 

In  this  report  we  take  advantage  of  the  fact  that  the  MPIE  is  expressed  in  terms  of 
potentials,  that  the  potentials  are  determined  by  a  gauge  condition  and  that  the  fields  are 
gauge  invariant  [2].  It  is  the  intention  of  this  investigation  to  show  that  integral  equation 
scattering  formulations  based  on  gauges  other  than  the  Lorentz  gauge  can  improve  numerical 
computation  speed  and  extend  the  range  of  applicability  of  contemporary  integral  equation 
scattering  models. 

In  classical  electromagnetic  field  analysis,  when  the  fields  are  expressed  in  terms  of  poten¬ 
tials,  the  Lorentz  gauge  historically  has  been  the  gauge  of  choice.  An  exception  to  this  rule 
is  sometimes  made  for  a  small  class  of  problems  where  the  divergence  of  both  the  magnetic 
vector  potential  A  and  the  current  density  J  are  zero  [3].  In  that  case  it  is  convenient  to 
use  what  is  commonly  referred  to  as  the  radiation  gauge  [2].  The  label  “radiation  gauge”  is 
relevent  l^ecause  the  conditions  described  above  respectively  include  and  infer  the  constraints 
V  •  A  =  0  and  $  =  0  that  are  associated  with  a  far  radiated  field. 

During  the  initial  stage  of  this  research  it  was  found  that  there  were  no  classical  Maxwell 
equation  scattering  models  based  on  the  Coulomb  gauge  ( Y  •  A  =  0,  $  ^  0),  except  for  one 
incorrect  attempt  described  in  a  text  by  Smythe  [4].  Because  the  Coulomb  gauge  offered  an 
alternative  to  the  Lorentz  gauge  and  because  it  was  as  yet,  an  unsolved  problem  in  scattering 
analysis,  it  was  decided  by  the  authors  that  a  thorough  investigation  of  the  Coulomb  gauge 
would  provide  the  necessary  “stepping  stone”  for  other  gauge  formulations.  Consequently 
the  remainder  of  this  report  will  detail  comparisons  between  the  Coulomb  and  the  Lorentz 
gauge  formulations  for  various  scattering  configurations  and  subject  to  a  variety  of  boundary 
conditions. 

A  change  in  gauge  essentually  means  that  the  vector  and  scalar  potentials  will  have 
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Green’s  functions  other  than  the  familiar  Lorentz  gauge  Green’s  functions.  Comparison 
between  the  Lorentz  gauge  and  other  gauges  can  best  be  displayed  by  incorporating  the 
effects  of  some  object  such  as  a  cylinder  or  wedge  in  the  Green’s  function,  thereby  leaving  the 
current  distibution  on  some  nearby  arbitrary  surface  as  the  only  unknown.  This  approach 
is  sometimes  referred  to  as  the  hybred  MM/Grecn’s  function  method  [5].  An  important 
consequence  of  the  hybred  MM/Green’s  function  method  presented  here  is  that  the  Green’s 
functions  are  in  their  most  general  (dyadic)  form. 

In  Chapter  2  magnetic  vector  and  scalar  potential  Green’s  functions  are  derived  subject 
to  the  Coulomb  gauge  for  a  general  directed  source  in  free  space.  It  is  shown  that  the 
Coulomb  gauge  vector  potential  Green’s  function  can  be  extracted  from  the  well  known 
free  space  Lorentz  gauge  vector  potential  Green’s  function  by  taking  advantage  of  a  set  of 
dyadic  identities.  This  method  avoids  the  difficulties  associated  with  deriving  dyadic  Green’s 
functions  from  first  principles  as  is  done  in  later  chapters. 

In  order  to  confirm  the  correctness  of  Coulomb  gauge  potential  expressions,  MPIE's  are 
derived  subject  to  both  the  Lorentz  and  Coulomb  gauges  for  a  finite  length  cylinder  scatterer. 
The  MPIE's  are  solved  for  the  cylinder  surface  current  by  the  method  of  moments.  Because 
the  current  is  unique  (i.e.  independent  of  gauge)  it  is  used  in  this  report  as  a  basis  for 
comparing  the  Coulomb  and  Lorentz  gauge  formulations  in  Chapters  2,  4  and  5. 

The  Coulomb  gauge  method  for  solving  source-excited  boundary  value  problems  in  closed 
(wa’.’oguide)  regions  is  presented  in  Chapter  -3.  Here  the  vector  eigenfunction  method  [6.7]  is 
used  to  obtain  the  dyadii..  Green’s  functions.  This  method  was  found  to  be  idealy  suited  for 
deriving  Coulomb  gauge  Green’s  functions  owing  to  the  particular  vector  nature  (solenoidal) 
of  the  Coulomb  vector  potential. 

In  Chapter  4  and  5  two  two-dimensional  open  region  scattering  problems  are  presented: 
plane  wave  scattering  from  an  arbitrary  shaped  conducting  scatterer  in  the  presence  of, 
respectively,  a  conducting  cylinder  and  a  conducting  wedge.  In  these  two  cases  the  Lorentz 
scalar  potential  is  expressed  as  an  infinite  series  whereas  the  Coulomb  scalar  potential  is 
closed  form.  A  closed  form  scalar  potential  is  seen  as  one  of  the  primary  advantages  of  the 
Coulomb  gauge  over  the  Lorentz  gauge  in  mixed  potential  formulations. 

Chapter  6  contains  our  conclusions  concerning  in  particular  the  Coulomb  gauge  as  com¬ 
pared  to  the  Lorentz  gauge  and  in  general  the  possibility  of  using  other  gauges  in  Maxwell 


equation  scattering  formulations.  Also  in  this  chapter  are  our  recommendations  for  future 
research  on  alternative  gauges. 
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Chapter  2 

Analysis  of  a  Wire  Scatterer 


In  this  chapter  the  magnetic  vector  and  scalar  potential  Green’s  functions  are  derived 
^  subject  to  the  Coulomb  gauge  for  a  general  directed  source  in  free  space.  Some  of  the  unique 

features  of  the  Coulomb  gauge  are  discussed  followed  by  a  presentation  of  the  MPIE’s  for  a 
thin,  straight,  finite-length  cylinder  scatterer  formulated  subject  to  the  Lorentz  and  Coulomb 
gauges.  The  chapter  concludes  with  a  comparison  of  Lorentz  and  Coulomb  gauge  current 
®  distribution  and  computation  time  results. 


2.1  Fields  in  Terms  of  Potentials 

The  time  harmonic  electric  and  magnetic  fields  are  related  to  the  magnetic  vector  poten¬ 
tial  A  and  scalar  potential  according  to 

E_  =  —ju;A  —  (2-1) 

H_  =  -  V  X  A  (2.2) 


with  the  time  convention  exp  {jujt)  suppressed.  In  terms  of  the  cylinder  surface  current 
density  J  and  the  dyadic  vector  and  scalar  Green’s  functions,  Q_  and  G<j)  respectively,  the 
magnetic  vector  and  scalar  potentials  are  [S]: 


A{l) 

=  J  G,  {l\l)  •  J-iL)  dv' 

v' 

(2.3) 

<I>(r) 

(2.4) 
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2.2  Lorentz  Gauge  Green’s  Functions  in  IVee  Space 

The  mathematical  form  taken  on  by  the  Green’s  functions  depends  upon  the  gauge  chosen 
when  the  differential  equations  for  the  potentials  are  derived.  For  example  the  Lorentz  gauge 


V  •  A  = 


(2.5) 


with  (2.1)  and  (2.2)  substituted  into  Maxwell’s  equation  V  x  H  =  J +  ia;e£,  results  in  a 
set  of  decoupled  potential  equations 


VM  +  =  -fiJ 


= 


V- J 

ja.'e 


where  Equation  (2.6)  motivates  the  dyadic  Green’s  function  equation 

=  -  ^(r  -  rO 


(2.6) 

(2.7) 


(2.8) 


where  ^  =  L5  {r  —  r').  I  =  xi  +  yy  +  Iz  is  the  identity  dyadic  and  (5(r  —  r')  is  the 
Dirac  delta  function.  Equation  (2.7)  effects  a  similar  scalar  Green’s  function  equation 


=  -  Sir  -  r') 

Equations  (2.8)  and  (2.9)  have  familiar  solutions  [8]: 
1  0  O' 

0  1  0  ^ =  LG{l\l') 


(2.9) 


Ga{l\l)  = 


G$(r|r') 


AttR 

Lo  0  1.' 

f-jkB. 


(2.10) 

(2.11) 


with  i?  =  [(x  -  x'Y  +  [y  -  y'f  +  (z  -  z'f]^  ■ 


2.3  Coulomb  Gauge  Green’s  Functions  in  Free  Space 


The  Coulomb  gauge 


V/1  =  0 


(2.12) 
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(2.13) 

(2.14) 


with  (2.1)  and  (2.2)  leads  to  the  set  of  Coulomb  potential  equations 
+  k^A  =  —  ^  =  —  fiJf 

JUJt 

Any  vector  J  can  be  expressed  solely  in  terms  of  its  solenoidal  (divergenceless)  J*  and 
lamellar  (zero  curl)  parts.  Mathematically  the  properties  of  these  two  parts  of  a  vector 
can  be  described  by  V  x  7^  0,  V  •  J"*  =  0,  V  x  =  0,  V  •  ^0.  The  superscript 

s  in  (2.13)  designates  a  solenoidal  current  density  Jf  that  results  when  the  lamellar  current 
j}  is  removed  from  the  total  current  density  J  =  ■  That  =  jcueV^  can  be 

shown  by  taking  the  divergence  of  (2.13)  and  then  applying  the  Coulomb  condition  (2.12). 
The  Green’s  function  equations  pertaining  to  the  Coulomb  potential  equations  (2.13)  and 

(2.14)  respectively  are 

V^e;^(r|r')  +  k^gfjr\r')  =  -  f  (r  -  r')  (2.15) 

V^G<j.(r|r')  =  -  Sit  -  r')  (2.16) 

Notice  that  because  the  d\’adic  delta  function  is  solenoidal,  the  Coulomb  Green’s  function 
must  also  be  solenoidal. 

Equation  (2.16)  is  simply  the  Poisson  equation  for  a  source  in  free  space.  The  solution 
to  (2.16)  is 

G*(Ei£')  =  ^  (2.17) 

The  dyadic  Green’s  function  method  for  solving  (2.15)  is  in  general  formatable  [6]  but 
in  this  particular  case,  where  the  source  resides  in  free  space,  we  can  make  some  observa¬ 
tions  that  will  lead  to  a  method  for  finding  the  Coulomb  gauge  Green’s  function  that  is 
understandable  to  those  readers  with  a  basic  knowledge  of  dyadic  vector  operations. 

First  we  observe  that  the  solenoidal  part  of  the  Lorentz  gauge  vector  potential  Green’s 
function  is  equal  to  the  Coulomb  gauge  vector  potential  Green’s  function.  This  can  be  shown 
by  first  decomposing  the  dyadics  on  the  left  and  right  hand  sides  of  (2.8)  into  solenoidal  and 
lamellar  parts  and  then  equating  just  the  solenoidal  components.  The  result  is  equation 

(2.15) . 

Since  the  Lorentz  gauge  Green’s  function  is  already  available  (equation  (2.10))  our  strat¬ 
egy  is  to  e.xtract  its  solenoidal  part.  This  can  be  done  by  first  substituting  the  identity 
—  VxVx^  into  a  second  identity  [9] 
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Mr  -  r'j  =  - 
“  Att  R 


(2.18) 


to  get 


S(r  —  r')  =  — 
=^-  ~  ’  47r 


-VV-  ^  UVxVx  ^ 


R 


R, 


(2.19) 


The  second  term  on  the  right  hand  side  of  (2.19)  is  the  solenoidal  delta  function  M.  Now 
^^(zilziO  c3.n  be  found  by  evaluating  [10] 

S;i(zili:')  =  I  -  r:')-a^{r;'\r;)dv"  (2.20) 


=  —  VxVx 

4x 


G(L''|r') 


,  /  0(L-|’l) 

\r  —  r" 


(2.21) 


where  the  r"  integration  extents  over  all  space.  After  some  tedious  integration  and  differen¬ 
tiation  we  find  that 


Sxx  S  xy  Sxz 

q  C  C 

^yx  ^yy  *^yz 

L  Szx  Szy 


^-jkR 

AttR 


=  g,G{r\r') 


(2.22) 


w 


.’ith 


{1-Q)L  -  (1-30 


RR] 

R? . 


Q 


1  -f  jkR  - 

(W 


jkR 


R  =  {x  -  x')£-  +  (y  -  y')y  +  (^  -  z')z 


(2.23) 

(2.24) 

(2.25) 


This  completes  our  presentation  of  the  Green’s  functions.  The  gauge  invarience  property 
of  the  electromagnetic  field  [2]  insures  us  that  the  fields  given  by  (2.1)  and  (2.2)  are  unique. 
The  equivalent  current  appearing  in  (2.3)  and  (2.4)  is  also  unique  (i.e.  independent  of  gauge) 
and  it  can  be  used  as  a  basis  for  comparing  these  two  gauges. 


2.4  Discussion 

One  of  the  apparent  delimmas  of  the  Coulomb  gauge  method  is  that  when  Green’s  method 
is  applied  to  equations  (2.13)  and  (2.15)  the  resulting  vector  potential  integral  contains  a 
solenoidal  current  J_‘ ,  rather  than  the  desired  surface  current  J.  We  can  immediately  resolve 
this  proldcm  by  replacing  the  solenoidal  current  with  the  surface  current; 
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(2.26) 


i4(r)  =  /z  I  £;^(r|r')  •  =  fi  j  a\{r\r^)  ■  l{rL)dv' 


This  can  be  done  because  solenoidal  and  lamellar  vectors  are  orthogonal  and  since  the 
Coulomb  Green’s  function  is  solenoidal,  the  integration  process  in  equation  (2.26)  will  auto¬ 
matically  remove  the  lamellar  part  of  the  surface  current.  We  emphasize  that  one  can  only 
numerically  solve  for  the  surface  current  because  the  solenoidal  and  lamellar  currents  can 
not  be  numerically  quantified  since  they  contain  non-physical  terms  that  exist  throughout 
space.  By  anology  this  last  point  is  illustrated  with  equation  (2.19)  where  the  two  terms  on 
the  right  hand  side  are  respectively  the  lamellar  and  solenoidal  delta  functions.  These  terms 
are  everywhere  non-zero  functions  of  position,  but  together  they  collapse  to  a  delta  function 
that  is  non-zero  only  at  r  =  r'. 

Some  contrasting  aspects  of  equations  (2.10)  and  (2.22)  clearly  indicate  that  the  Lorentz 
gauge  has  a  computational  advantage  over  the  Coulomb  gauge.  For  example  the  Coulomb 
gauge  Green’s  function  contains  a  fully  populated  dyadic  coefficient  matrix  ^  while  the 
Lorentz  gauge  coefficient  is  simply  the  identity  dyadic  L  This  means  for  example  that  current 
flowing  in  the  z-direction  will  produce  three  (x,  y  and  z)  dissimilar  Coulomb  vector  potential 
components  or  a  single  (r)  Lorentz  vector  potential  component.  A  general  bent  wire  scatterer 
solved  b}-  the  method  of  moments  MPIE  technique  would,  for  most  of  the  moment  matrix, 
require  three  calculations  (a  column  of  the  ^  dyadic)  per  matrix  element  with  the  Coulomb 
gauge  and  only  one  calculation  with  the  Lorentz  gauge. 

On  the  other  hand  for  a  straight  wire  a  Coulomb  gauge  analysis  appears  at  first  glance 
to  have  some  advantages  over  the  Lorentz  gauge  because  of  the  simple  mathematical  form  of 
the  Coulomb  scalar  potential.  Although  this  will  not  prove  to  be  the  case,  we  present  below 
an  MPIE  analysis  of  a  straight  finite  length  cylinder  scatterer  illuminated  by  a  plane  wave 
because  it  is  a  simple  instructive  problem  for  comparing  our  gauge  formulations. 


2.5  Results  for  Wire  Scatterer 

A  straight  cylinder  of  radius  a  <<  A  and  length  £  is  illuminated  by  a  z-polarized  plane 
wave  traveling  in  a  direction  normal  to  the  cylinder  (z-)  axis.  By  enforcing  the  condition 
that  the  tangential  component  of  the  total  electric  field  is  zero  on  the  cylinder  surface  one 
gets  the  MPIE 
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where  in  the  Lorentz  and  Coulomb  gauges  the  potentials  are 


>1. 


G{z,z',ip')  1  ,  /  ,  /  /  Lorentz  1 

Coulomb} 


$  = 


27ra;e 


J>r 


G{z,z',(f')' 


1 

4-rR 


dip'  dz' 


Lorentz 

Coulomb 


(2.28) 


(2.29) 


R  =  [{z  —  z'Y  -\-2a^  sin^  9']^  and  represents  the  incident  electric  field.  Because  it  is  well 
documented  [1]  we  will  not  elaborate  here  on  the  method  of  moments  numerical  technique 
for  computing  the  current  on  a  thin  cylinder  scatterer  except  to  say  that  in  our  procedure 
we  employed  20  rectangular  basis  and  testing  functions. 

Fig.  2.1  shows  the  current  distributions  on  three-quarter  and  full  wavelength  thin  cylin¬ 
ders.  The  solid  lines  depict  real  and  imaginary  parts  of  the  current  due  to  a  Coulomb  gauge 
formulation  while  the  overlying  triangles  are  from  the  Lorentz  gauge  analysis.  The  two  re¬ 
sults  are  numerically  indistinguishable.  The  computation  time  difference  for  each  of  these 
two  results  was  15%  in  favor  of  the  Lorentz  gauge. 

The  far  zone  electric  field  can  be  found  by  letting  r  — »  00  in  (2.28)  and  (2.29)  and  then 
substituting  the  result  into  equation  (2.1).  The  far  electric  field  is 

Eeir)  =  j  (2.30) 

V 

in  both  gauges.  Because  the  currents  and  the  far  field  expressions  are  the  same,  there  is  no 
difference  in  the  far  fields  obtained  by  the  two  gauge  formulations. 

The  difference  in  the  two  gauges  emerges  in  the  vector  and  scalar  potential  p'  integrations. 
The  dyadic  Szz  is  a  regular  non-singular  function  that  presents  no  numerical  difficulties.  The 
singularity  in  G{z^z' is  evaluated  by  substracting  and  adding  1/R  as  follows: 


J  G  {z,z',p')dp'  = 


-  1  ,  ,  1 


(2.31) 


so'  =  i/>'=:-7r 

Evaluation  of  the  Lorentz  scalar  potential  requires  both  integrals  on  the  right  hand  side 
(RHS)  of  (2.31)  while  only  the  second  RHS  integral  is  required  for  the  Coulomb  scalar 


potential.  This  advantage  gained  with  the  Coulomb  gauge  in  the  scalar  potential  calculation 
is  not  enough  to  overcome  the  fact  that  the  integration  in  the  Lorentz  gauge  can  be  done 
once  and  then  used  in  both  the  scalar  and  vector  potentials  (c.f.  equations  (2.28)  and  (2.29)). 
We  can  then  conclude  that  the  Lorentz  gauge  potentials  can  be  computed  more  rapidly  than 
the  Coulomb  gauge  potentials  for  all  thin  wire  configurations  and,  it  is  safe  to  say,  in  any 
case  where  a  wire  is  in  the  presence  of  an  object  for  which  a  Green’s  function  can  be  derived. 


1.5 


Figure  2.1;  Comparison  of  Lorentz  and  Coulomb  gauge  current  distribution  on  a  (c 
-  dipole  and  a  (b)  A  -  dipole.  The  incident  electric  field  is  set  to  unit}'. 


)  A/4 
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•  Chapter  3 

Source  in  a  Waveguide 


There  are  two  general  categories  for  electromagnetic  boundary  value  problems:  wave  scat¬ 
tering  and  wave  guidance.  In  this  chapter  we  present  a  Coulomb  gauge  analysis  of  the  fields 
in  a  source-excited  waveguide.  In  our  waveguide  example  the  vector  eigenfunction  expansion 
method  is  used  which  is  essentially  the  same  as,  but  mathematically  less  complicated  than, 
the  wave  scattering  analysis  found  in  succeeding  chapters. 

Historically  it  was  thought  that  the  scalar  and  vector  potentials  were  separable  only  in  the 
Lorentz  gauge.  However  it  is  shown  in  this  chapter  that  for  perfectly  conducting  boundaries 
the  scalar  and  vector  potentials  are  separable  in  the  Coulomb  gauge  by  taking  advantage 
of  the  solenoidal  and  lamellar  properties  of  vectors  as  was  done  in  principle  for  a  source  in 
infinite  free  space  in  Chapter  2. 


3.1  Coulomb  Gauge  Vector  Potential  Green’s  Func¬ 
tion  for  a  Rectangular  Waveguide 


Consider  a  perfectly  conducting  rectangular  waveguide  aligned  along  the  2-axis  and  with 
dimensions  a  and  b  along  the  x-  and  j/-axes,  respectively.  To  find  subject  to  the  boundary 
condition  that  n  x  =  0  on  the  waveguide  interior  surface  we  follow  the  Ohm-Rayleigh 
method,  as  described  by  Tai  [6].  First,  we  expand  in  (2.15)  as  follows: 


^  oo  oo 

d’(Llr')  =  dhZT. 


m  =  0  n=0 


2-^0 

Trabk^ 
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•  (3.1) 

/c 

where  and  are,  respectively,  the  even  and  odd  solenoidal  eigenfunctions  of  the 

operator  V  x  Vx,  satisfying  the  Dirichlet  boundary  conditions  at  the  waveguide  walls. 
These  vector  eigenfunctions,  which  share  the  eigenvalues  /c^,  are  generated  by  the  scalar 
eigenfunctions  of  t^e  operator  —  according  to  [6] 

(3.2) 

=  YxM^^mnih)  (3.3) 

where  ipemn  and  ipomn  satisfy,  respectively,  the  Neumann  and  Dirichlet  boundary  conditions 
at  the  waveguide  walls.  For  the  rectangular  waveguide,  one  easily  finds 


cos  ^  cos  ^ 
sin  ^  sin  ^ 

a  0 


4’‘mn{h)  = 


+  e  =  Ht) 

The  primed  functions  M'  and  iV'  in  (3.1)  depend  on  the  source  coordinates  x',y',  and  z'. 
The  meaning  of  the  symbol  So  in  (3.1)  is  as  follows:  <5o  =  1  when  m  =  0  or  n  =  0  and  =  0 
otherwise. 


With  (3.1)  in  mind,  we  now  expand  as  [6] 

°°  9  _  X 


m=0  n=0 


■  +  ■  (3.7) 

K  J 

To  evaluate  the  coefficients  a{h)  and  b{h)  we  substitute  the  expansions  (3.1)  and  (3.7)  into 
(2.15),  introduce  the  operator  under  the  integration  and  summation  signs  on  the  left  side 
and  make  use  of  the  relation  and  a  similar  relation  for  the  TV  functions.  As 

a  result,  we  find 


a{h)  =  b{}i)  = 


/i2  - 


where 


—  k],  kc  <  k 

kc>k. 
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We  can  now  express  (3.7)  in  the  form 


aAir-lr!)  = 


1  <50 


U2„}y  ^  _P 

“^m=0n=0 


°°  flh 


+ 


dh .  kl  kl 


f  ^ 
J  h?^K^- 


+ 


k]  h^^kl 


)ILUh)N!^A-h) 


(3.10) 


which  clearly  exhibits  two  sets  of  poles  in  the  complex  /i-plane,  located  at  ±.kg  and  at  ±j/:c 
(h  =  0  is  a  regular  point  of  the  integrand).  The  integrals  in  (3.10)  can  be  easily  evaluated 
by  residue  calculus  and  the  result  can  be  expressed  as 

fi^(rlr')  =  ilC,(r|r')  -  fi,(r|r'))  (3.11) 

where 


a,{L\r:)  = 
+ 


1  oo  oo  _  c 

n  E  E  -^lk^iL„.i±K)M:,u^k,) 

m=0n=0 

iL^n{±k,)£CU^k,)l  z<z' 


(3.12) 

(3.13) 


and 


^  E  E  ziz'. 

m=l  n=  1 


We  observe  that  is  the  static  limit  of  i.e., 


(3.14) 


Ocitlt)  =  lim£j(r|r').  (3.15) 

Finally,  we  remark  that  the  integral  in  (3.1)  can  also  be  evaluated,  with  the  result 
i’iz -  l)  =  (ii  +  iy)<5(!: -  l)  +  ^(r|r').  (3.16) 

3.2  Coulomb  Gauge  Scalar  Potential  Green’s  Func¬ 
tion  for  a  Rectangular  Waveguide 

The  dyadic  Green’s  function  for  the  lamellar  part  of  the  electric  field  can  be  assembled 
in  a  few  simple  steps  beginning  with  the  solution  of  the  static  problem 
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(3.17) 


V^G^(rfc')  =  -%-!:') 

where  must  vanish  at  the  waveguide  walls  and  at  infinity. 


For  the  rectangular  waveguide,  we  easily  find  [4] 

A 

ab 


ry  /  \  2  ~  ~  1  .  rmrx'  .  nny'  .  mnx  .  mry  _kAz~z'\ 

G^{r\r)  =  -  ^-sin-— sin—  sm  — sin— e  ' 


where  is  defined  in  (3.5). 

Finally,  we  observe  that 

-  r')  =  1 1  <5(r  -  r')  -  ^(dn)- 


3.3  Discussion 


(3.18) 


(3.19) 


The  attractive  feature  of  the  Coulomb  gauge  is  the  explicit  separation  of  the  electric  field 
into  its  lamellar  and  solenoidal  constituents.  The  lamellar  part,  which  contains  the  dominant 
R~^  singularity,  where  R  is  the  distance  between  the  source  and  the  observation  points,  is 
easy  to  determine,  since  the  scalar  potential  $  can  be  obtained  by  a  simple  differentiation  of 
the  corresponding  static  potential.  (For  a  few  simple  geometries,  this  potential  can  be  ob¬ 
tained  in  closed  form  by  image  theory).  The  remaining  part  comprises  the  more  manageable 
R~^  singularity  and  is.  of  course,  solenoidal  everywhere,  including  the  source  region.  Hence, 
if  the  eigenfunction  expansion  technique  is  employed,  jE  can  be  conveniently  represented  in 
terms  of  only  the  solenoidal  M  and  iV  functions,  and  the  lamellar  L  functions  are  indeed 
obviated  [6]. 

The  price  paid  for  these  advantages  of  the  Coulomb  gauge  is  the  added  difficulty  in 
solving  for  the  vector  potential  if  the  approach  is  taken  of  expanding  A  (or  E_  in  terms  of  the 
E  and  H  waveguide  modes  [3,  11,  12].  This  difficulty  is  due  to  the  fact  that  (2.13),  unlike 
the  corresponding  equation  in  the  Lorentz  gauge,  involves  the  solenoidal  part  of  J  which  is 
usually  a  much  more  complicated  function  than  J  itself.  For  example,  J  associated  with  a 
point  dipole  has  the  simple  form  of  the  Dirac  delta,  whereas  the  corresponding  Jf  and  Jf  are 
not  localized  at  a  single  point  in  space.  These  difficulties  in  obtaining  the  Coulomb  vector 
potential  can  perhaps  be  blamed  for  the  subtle  error  in  Smythe’s  book  [4],  which  was  written 
years  before  the  nature  of  the  field  in  the  source  region  was  fully  explored. 
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Chapter  4 


Transverse-Electric  Wave  Scattering 
by  Two-Dimensional  Surfaces  of 
Arbitrary  Shape  in  the  Presence  of  a 
Circular  Cylinder 


In  this  chapter  we  return  to  the  open  region  problem  but  here  we  construct  a  hybrid 
MM/Green‘s  function  problem.  Two  forms  of  the  MPIE  are  developed  for  two-dimensional 
perfectly  conducting  (PC)  surfaces  of  arbitrary  shape  in  the  presence  of  an  infinite  PC  cylin¬ 
der  of  circular  cross-section,  subject  to  transverse-electric  excitation.  One  of  the  MPIEs  is 
based  on  the  Coulomb  gauge;  the  other  imploys  the  Lorentz  gauge.  In  either  case,  the  effect 
of  the  cylinder  is  incorporated  in  the  integral  equation  by  means  of  the  appropriate  Green’s 
functions,  leaving  the  current  distribution  on  the  arbitrary  surface  as  the  only  unknown. 
The  Green’s  fuucticns  are  derived  by  the  eigenfunction  expansion  technique.  An  existing 
well-established  moment  method  procedure  is  adapted  to  numerically  solve  both  forms  of 
the  MPIE.  Computed  results  are  presented  for  several  cases. 


4.1  Statement  of  the  Problem  and  Assumptions 


The  geometry  of  the  problem  considered  is  shown  in  a  cross-sectional  view  in  Fig.  4.1, 
which  depicts  a  surface  S  of  arbitrary  shape  (open  or  closed)  near  a  circular  cylinder  of 
radius  a,  whose  axis  coincides  with  the  z  axis  of  a  Cartesian  coordinate  system  {x,y,z)  with 
an  inscribed  cyhndrical  system  (p,^,  2).  The  infinitesimally  thin  surface  and  the  cylinder 
are  made  of  perfect  conductors  and  are  of  infinite  extent  and  invariant  along  the  2  coordi¬ 
nate.  The  surrounding  medium  is  homogeneous  and  is  characterized  by  permittivity  e  and 
permeability  /z.  The  structure  is  illuminated  by  a  plane  wave  transverse  electric  (TE)  to 
the  cylinder  axis.  .An  time  convention  is  assumed  and  suppressed.  The  quantities  of 
interest  are  the  current  distribution  on  the  arbitrarily-shaped  surface,  the  far  zone  field,  and 
the  bi.static  radar  cross-section  (RCS). 


4.2  Lorentz  and  Coulomb  Gauge  Green’s  Functions 
for  the  Cylinder 


do  determine  we  employ  the  eigenfunction  expansion  technique  [6],  in  which  the 
solution  is  sought  in  terms  of  the  eigenfunctions  of  the  vector  Laplacian.  These  vector 
eigen  1  unctions  can  in  turn  be  expressed  in  terms  of  two  sets  of  scalar  eigenfunctions  a-nd 
which  satisfy  the  equaticm 


(4.1) 


where  V/  is  the  transverse  Laplacian  in  cylindrical  coordinates,  subject  to,  respectively,  the 
Dirichlet  and  the  Neumann  boundary  conditions  at  p  =  a.  Other  conditions  that  must  also 
be  imposed  on  V’nrt  and  are  the  periodicity  in  (p  and  boundedness  as  p  — >  00.  To  find 
V’nK  and  the  eigenvalues  we  follow  the  separation-of-variables  procedure,  which  yields 


r„(Ka) 

fn(«o) 


(4.2) 


wheif'  ri  =  0,±1  ±  2,...,  and  0  <  k  <  oc.  In  the  last  equation,  and  denote, 
respi'i  t ivrly,  the  Bessel  function  and  the  Hankel  function  of  the  second  kind,  both  of  order 
n.  and  L,,  and  Fn  denote  the  "reflection  coefficients,’’  given  as 


(4.3) 


r„(Ka)  -  ,2),  r„(«:a)  -  ,^y 

Hn  («a)  Hn  {ku) 

where  and  denote  derivatives  with  respect  to  the  arguments  of  the  respective  func¬ 
tions.  It  can  be  shown  that  t/)„«  and  ‘ipnK  possess  the  orthogonality  properties 

The  set  of  the  vector  eigenfunctions  consists  of  three  species  [6,  13],  M.nni  ‘‘•nd  N_nK^ 
where 


Llixk  —  X-'-'nx 


Al-yiK  3Z  ^  ±.'^nK  ^  ^  y^riH 


'X.nw  =  ■■  X  V  X  Zxl^riK  —  Ixi/-’,, 


Clearly,  the  functions  are  lamellar  and  the  and  functions  are  solenoidal  and 
it  can  be  shown  that  all  of  them  satisfy  (4.1).  From  (4.4)  through  (4.7),  it  follows  that  the 
vector  eigenfunctions  possess  the  orthogonality  properties 

.  „  f  u  1  f  i;.<.  1 


J  J  j  ALnK  f-j  Kn'„'  }pdpdip  =  Snn'^c'S{K-H')  (4.S) 

We  now  expand  in  (2.8)  in  terms  of  the  vector  eigenfunctions  and  in%'oke  (4.1) 
and  (4.8)  to  obtain 


^  OC’  , 

r  fj  pr 

fl'4(£i£'i = j:J 


oc  ^  / 

a(^ie'i=  s:  / [M..(g)M;.(£')  +  iv..(g)ivt.(g': 


(4.10) 


Noting  that  the  last  term  in  (4.10)  contributes  only  the  zz  component  of  G'j  (which  is  not 
relevant  to  tin-  IE  problem  considered  here),  we  separate  it  by  writing 
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(4.11) 


a\  =  a\,+iiG‘, 


where  comprises  only  components  transverse  to  z. 
Upon  substituting  (4.9)  into 

2-S^(£l£')  =  -V'G»(£l£') 


(4.12) 


and  using  (4.5)  and  (4.1),  we  readily  obtain  the  Lorentz  gauge  scalar  potential  Green’s 


function  as 


Gitip  \i)=  £  7 (p)  Ck  {pl) 


(4.13) 


In  view  of  (4.7),  it  immediately  follows  that  =  6'$. 

The  integrals  encountered  in  (4.9),  (4.10),  and  (4.13)  can  be  evaluated  by  the  residue 
calculus  [13,  14].  Hence,  we  obtain 


OC 

J  R„(kp)R-Jkp')  — 


(»!  -  k'^) 


n-  [ 


(4.14) 


=  ^{li„(.£n//-(l-£>)-;^|(g)”-g)"]}  (4.15) 

and 

. 

/  RnMR:M^^^  =  (4.16) 

0 

where  p<  and  p>  denote,  respectively,  the  lesser  and  the  greater  of  p  and  p'.  The  first  right 
side  member  in  each  of  these  equations  is  contributed  by  the  residue  of  a  pole  at  k  =  k 
of  the  respective  integrand.  The  remaining  right  side  terms  in  (4.14)-(4.15)  arise  from  the 
integration  around  the  singularities  at  the  origin  of  the  K-plane.  These  “static”  terms  have 
been  overlooked  in  some  past  work,  as  pointed  out  in  [13]. 

When  (4.14)-(4.16)  are  substituted  into  (4.9)-(4.10)  and  (4.13),  we  expectedly  find  that 
the  resulting  expressions  can  in  each  case  be  broken  into  two  series,  of  which  one  is  ana¬ 
lytically  summable  and  represents  the  “primary”  (free-space)  part  of  the  respective  Green’s 


24 


function,  and  the  other,  which  must  be  summed  numerically,  accounts  for  the  presence  of 
the  cylinder.  We  will  distinguish  these  two  parts  by  superscripts  “p”  and  “c,”  respectively. 
Hence,  with  the  notation 


R=\p-  p'\^  i  = 


Bi  =  Hl^\kR)  - 
B2  =  h^2\^R)  - 


-K{kRy 


we  can  express  the  primary  components  of  as 


=  L  _  ^2  — sin^^ 

'  kR  R?  ^ 


=  ij  -  P'cos()s,na 

G^.  =  i  B.  ^  +  A  i(/)'  -  p  cos  f  )(/>  -  P'  cos  f ) 

The  remaining  secondary  components  take  the  form 

GZ'  =  r  E^cosnf  lj^t„(kp)n!.»(kp'}/i!,^>(kp)  -  oJ 

^  [k^pp' 

1  00  j-f 

G;Z  =  ^  E  2sin  ’<(  ^t„(ka)//<nkp')N<‘l(kp)  +  (?, 

‘>S  I*-’'’  J 

1  00  r  ^ 

G’Z’  =  -^^2smn(  ;^r„(io)Bf  (M  +  Qn 

4;  [A-r/ 


=  P  £  e„cos7i^  [r„(A:a)//f’  (M  “ 

4./  n=0 

/here  e„  =  1  for  77  =  0  and  £„  =  2  for  n  >  0  (Neumann’s  number),  and 


Qn  = 


]n  /  a 


-[kaY  \pp' 


2  \  n  +  l 


(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 

(4.26) 

(4.27) 

(4.28) 


This  last  factor,  as  well  as  the  second  right  side  members  in  (4.18)-(4.19),  are  the  previously 
mentioned  static  terms. 

Finally,  the  primary  and  secondary  parts  of  the  Coulomb  scalar  potential  Green’s  function 
Go  are  found  in  closed  form  as  ([15],  p.  69) 

Gl  =  ~^\nR  (4.29) 

Ztt 

(4.30) 

Observe  that,  although  we  are  solving  a  time-harmonic  problem,  this  Green’s  function  is 
“static,”  he.,  it  does  not  incorporate  the  time  retardation  effect.  This  paradox  is  character¬ 
istic  of  the  Coulomb  gauge  [16]. 

We  can  obtain  the  lamellar  part  of  the  vector  potential  Green’s  function  by  a  procedure 
analogous  to  that  followed  above  for  the  solenoidal  part.  Adding  the  so-obtained  to 
the  previously  found  results  in  the  Lorentz  vector  potential  Green’s  function  which 
we  also  decompose  into  the  sum  of  the  primary  part  and  the  secondary  part  The 


transverse  components  of  the  former  are  found  as 

G"/  =  =  ^H!i‘HkR)cos(  (4.31) 

G^„.  =  -G%,.  =  jrH^^\kR)  sin  {  (4.32) 

while  the  corresponding  components  of  the  latter  are  given  as 

1  OO 

=  3-  S  \r„(k„)Hj.^r(ki,')H!^y(kp) 

n=0 

^-^r,(ka)H'-:\kp')H^(kp)  (4.33) 

1  OO  r  „ 

GV  =  P  E  7-;r,(ia)//ra(t/)Hf''(i-/>) 

n=l 

+  Jr.(ta)«W(V)//t">(M  (4.34) 

1  OO  r 

G't.'  =  -  E  h^r,(t«)//f '(V)//'">(M 

+  ~t,(ka)Hy^\kp')H'^y{kp^  (4.3.5) 
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1  °° 
n=0 


n* 


k^pp 


■>rUka)H^''\kp‘)H^^y(kp) 


(4.36) 


It  turns  out  that  the  static  terms  of  Qf^  and  cancel  when  these  dyadics  are  added  to 
form  which  was  previously  observed  in  [17]. 

Upon  substituting  (4.16)  into  (4.13),  we  also  readily  find  the  primary  and  secondary 
parts  of  the  Lorentz  scalar  potential  Green’s  function  as 


Gl  = 

V 

1  oo 

enCosn(Tr,{ka)Hl^\kp')H^^\kp) 

71  =  0 


4.3  Comments 


(4.37) 

(4.38) 


We  observe  that  the  primary  terms  (4.3l)-(4.32)  of  the  Lorentz  gauge  Green’s  function  are 
much  simpler  than  the  corresponding  terms  (4.20)-(4.23)  of  its  Coulomb  gauge  counterpart. 
.^Iso,  which  is  perhaps  not  obvious,  the  diagonal  {pp'  and  (pp')  primary  components  of  the 
latter  exhibit  a  logarithmic  singularity  as  R  0,  whereas  the  off  diagonal  (pp'  and  pp') 
components  are  well-behaved  there.  (The  same  is  clearly  true  about  the  Lorentz  gauge  vector 
potential  Green’s  function).  We  note  that  were  it  not  for  the  static  terms  (the  second  right 
side  members  in  (4.18)-(4.19)),  and  would  have  nonintegrable  singularities  at  the 
source  point. 

The  static  terms  Qn  also  help  the  convergence  of  the  series  (4.24)-(4.27),  which  would 
otherwise  diverge.  However,  even  with  this  help,  the  large-n  behavior  of  the  terms  of  these 
series  is  0(n°).  In  contrast,  the  terms  of  the  corresponding  Lorentz  gauge  series  (4.33)-(4.36) 
behave  as  0(n~^ )  for  large  n — which  is  better,  but  not  good  enough.  Consequently,  the  series 
of  both  formulations  were  accelerated  by  Kummer’s  transformation  ([18].  p.  16),  resulting 
in  series  with  terms  behaving  as  0(n~^).  Because  this  procedure  has  now  become  standard 
([10],  p.  819,  [19],  [20]),  we  omit  the  details  here. 

Finally,  we  remark  that  there  are  two  procedures  that  one  could  follow  to  derive  the 
dyadic  Green’s  functions  for  the  problem  of  Fig.  4.1  —  the  so-called  scattering  superposition 
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method  described  in  [6]  and  the  method  of  direct  eigenfunctions  [17],  which  from  the  outset 
incorporates  the  effect  of  the  cylinder.  We  have  chosen  to  present  the  latter  approach  here, 
because  we  had  difficulty  justifying  the  scattering  superposition  technique  when  the  static 
terms  are  present. 

4.4  Numerical  Method 

To  numerically  solve  the  MPIEs  in  both  gauges,  we  adapt  the  MM  procedure  developed 
in  [21]  for  objects  in  free  space.  In  this  approach,  the  arbitrary  surface  is  modeled  in  terms 
of  straight  line  segments,  and  piece-wise  linear  expansion  functions  (he.,  triangle  functions 
which  span  two  adjacent  segments)  are  employed  to  approximate  the  current.  The  same 
functions  are  used  to  “test”  the  integral  equation,  producing  an  algebraic  system  of  equations, 
which  is  readily  solved  for  the  current  expansion  coefficients.  The  procedure  given  in  [21] 
had  to  be  properly  modified  to  account  for  the  dyadic  nature  of  the  vector  potential  Green’s 
functions.  We  omit  the  details  in  the  interest  of  brevity. 

4.5  Scattered  Far  Fields 

The  “incident”  field  due  to  a  TE  plane  wave  with  the  amplitude  can  be  easily  found 
([22],  p.  2-35)  and  is  not  listed  here  for  the  sake  of  brevity. 

The  far  zone  magnetic  field  can  be  found  by  using  either  the  Lorentz  or  the  Coulomb 
vector  potential  in  (2.2)  and  letting  p  —y  oo.  The  result  is 

(4.39) 

with  Hq  =  Eq/t],  where  ?;  is  the  intrinsic  impedance  of  the  medium. 

The  remaining  terms  in  (4.39)  are  defined  as 

E^ip',0  =  (4.40) 


(4.40) 


2S 


(4.41) 


£'(/>'.  0  =  (V)sinnf 


n=l 


£  ^nr'''^t'n{ka)H^^^'{kp')  cos  ni 


n=0 


^(0  =  E^"(-l)"rn(A:a)cosn^ 


(4.42) 


n=0 


The  first  and  the  second  terms  in  the  curly  brackets  in  (4.39)  represent,  respectively,  the 
“direct”  and  the  “reflected”  (by  the  cylinder)  parts  of  the  far  field  produced  by  the  current 
J  on  5,  while  the  last  term  is  the  far  zone  “reflected”  field  produced  by  the  incident  wave. 
Finally,  the  bistatic  RCS  a  is  defined  as  ([22],  p.  358) 


<T  ,  I 


(4.43) 


and  is  only  a  function  of  ip  and  y?,. 


4.6  Validation  and  Sample  Results 

We  present  in  this  section  sample  numerical  results  for  the  structures  illustrated  in 
Fig.  4.2.  These  simple  structures  have  been  chosen  for  analysis  to  validate  the  computer 
programs  based  on  the  Lorentz  and  Coulomb  gauges  and  to  compare  their  efficiency. 

In  all  four  cases  shown  in  Fig.  4.2  the  cylinder  radius  a  =  0.25A  {ka  —  7r/2),  with  A  =  1  m. 
In  the  problem  of  Fig.  4.2(a)  the  cylinder  resides  at  the  center  of  an  equilateral  triangular 
box  and  is  completely  shielded  from  the  incident  field.  Therefore,  its  presence  is  irrelevant  to 
the  exterior  scattering  problem,  and  we  may  compare  our  results  with  that  obtained  for  the 
empty  triangular  box  by  a  free  space  code  [21].  Thus,  in  Figs.  4.3(a)  and  4.3(b)  are  shown, 
respectively,  the  current  distribution  on  the  triangular  box  and  the  bistatic  RCS  (4.43)  for 
pi  =  90°,  obtained  by  both  methods.  Good  agreement  between  the  corresponding  results 
is  observed.  Because  the  Coulomb  and  Lorentz  gauge  based  programs  produced  graphically 
indistinguishable  results,  only  the  Lorentz  gauge  data  are  displayed  for  this  and  other  case 
discussed  here. 

It  is  of  interest  to  note  that  for  the  geometry  of  Fig.  4.2(a)  the  last  two  terms  of  (4.39) 
must  cancel  each  other,  because  the  far  field  is  completely  specified  by  the  first  right  side 
term.  When  the  computed  approximate  J  is  used  in  (4.39),  this  cancellation  is  not  perfect, 


which  is  the  likely  reason  behind  the  slight  discrepancy  in  the  RCS  results  observed  in 
Fig.  4.3(b). 

The  remaining  structures,  illustrated  in  Figs.  4.2(b)  through  4.2(d),  are  variations  on 
the  first  problem.  The  corresponding  results  for  two  angles  of  incidence,  </?,■  =  60°  and  90°, 
shown  in  Figs.  4.4  through  4.6,  illustrate  the  effect  of  removing  parts  of  the  triangular  shield 
on  the  current  distribution  and  the  RCS.  In  these  examples,  a  15%  savings  in  the  matrix  fill 
time  was  realized  when  the  Coulomb  gauge  was  used  instead  of  the  Lorentz  gauge. 
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Figure  4.2:  .\  test  problem  in  which  the  surface  5  of  Fig.  4.1  forms  a  triangular  box  end' >?in; 
the  cylinder,  (b)-(d)  \'ariations  of  (a)  in  which  part.^i  of  the  triangular  shield  are  removed 
In  all  cases  A  =  Im. 
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Figure  4.5:  (a)  Normalized  current  density  on  5  and  (b)  bistatic  RCS  for  the  configuration 
of  Fig.  4.2fc)  with  tp,  =  60°  and  90°. 
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r'iguif'  4.6: 


(a)  Normalized  current  den.sily  on  S  and  (b)  bistatic  RCS  for  the 


configuration 


N  Fic.  4.2(d)  with  =  60°  and  90°. 


Chapter  5 


Transverse-Electric  Wave  Scattering 
by  Two-Dimensional  Surfaces  of 
Arbitrary  Shape  in  the  Presence  of  a 
Wedge 


Here  we  investigate  an  infinitesimally  thin,  perfectly  conducting  (PC)  surface,  which  is 
of  infinite  extent  and  invariant  along  the  z  axis,  but  has  an  arbitrary  (open  or  closed)  cross- 
sectional  profile  C,  and  which  resides  near  the  edge  of  an  infinite,  PC  wedge  of  angle  o ,  whose 
upper  surface  coincides  with  the  (a;  >  0,  y  =  0)  half-plane  (see  Fig.  5.1).  The  surrounding 
medium  is  homogeneous  and  is  characterized  by  permittivity  e  and  permeability  y.  The 
structure  is  illuminated  by  a  plane  wave  transverse  electric  (TE)  to  the  z  axis.  We  consider 
scattering  from  a  wedge  because  it  is  one  of  the  canonical  problems  of  classical  field  theory. 

5.1  The  Wedge  MPIE 

In  the  case  considered  here  (Fig.  5.1 ),  J  is  the  unknown  current  induced  on  the  PC  surface 
C  by  a  known  ‘‘incident”  field  E’,  which  exists  in  the  presence  of  the  wedge.  Enforcing 
the  condition  that  the  tangential  components  of  E  and  E*  cancel  on  C  leads  to  the  EFIE 


where  r  is  a  unit  vector  normal  to  C  at  £.  This  equation  is  then  solved  for  J_  by  the  MM 
procedure. 

5.2  Lorentz  and  Coulomb  Gauge  Green’s  Functions 
for  a  Wedge 

We  begin  by  constructing  two  scalar  eigenfunctions,  0n«  and  ipuK,  which  satisfy  the  equa- 


^nx{p) 

'PnK{p) 


subject  to,  respectively,  the  Dirichlet  and  Neumann  boundary  conditions  on  the  PC  wedge. 
They  must  also  satisfy  the  edge  condition  and  be  bounded  as  p  — >  oo.  (In  (5.2),  is  the 
transverse  Laplacian  in  cylindrical  coordinates  (no  z  dependence).)  These  eigenfunctions  are 
easily  obtained  by,  for  example,  the  characteristic  Green’s  function  procedure  [23].  For  the 
problem  of  Fig.  5.2,  the  normalized  eigenfunctions  are 


Pn^ip) 

PnK{p) 


P  I  COS  l/n'-pf 


where  is  the  Neumann  number  (e„  =  1  for  n  =  0  and  €„  =  2  for  n  >  0),  is  the  Bessel 
function  of  fractional  order  Vn  =  nirll3^  y5  =  27r  — a,  0<«<  oo,  and  n  =  0,1,2,...  These 
eigenfunctions  possess  the  orthogonality  properties 

where  8nn'  is  the  Kronecker  delta  ((5„„/  =  1  when  n  =  n'  and  zero  otherwise). 

We  now  construct  three  species  of  vector  eigenfunctions,  M.nKi  K-tiki  ^-s  [6,10] 

Ln^  =  (5-5) 

M.r,«  =  X  i0n«  =  -1  X  (5.6) 


KrxK  =  -5.  X  V  X  =  iKV’n. 

K 
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Clearly,  the  L^k  functions  are  lamellar  and  the  and  functions  are  solenoidal.  It  is 
easily  confirmed  that  all  of  them  satisfy  (5.2)  and  the  required  boundary  conditions.  It  also 
follows  from  (5.4)  that  the  vector  eigenfunctions  possess  the  orthogonality  properties 


P  oo 

// 

0  0 


L  1 

{  L  ,  ,  \ 

Hn'K'  1 

> .  < 

.  N^'.>  J 

(5.8) 


VVe  now  expand  in  (2.8)  in  terms  of  the  vector  eigenfunctions  and  employ  (5.2)  and 

(5.8)  to  obtain 


_  OO  , 

^  dK 


2!.  (£  I  £')  =  Z  / 


(5.9) 


and 


^  OO  . 

/*  (Ik.  r 

I  e')  =  E  /  [m„.(£)m„,(p')  +  w„.(£)iv„.(£') 


n=0  ) 


c(/c2  -  P) 


(5.10) 


Noting  that  the  term  in  (5.10)  comprising  the  functions  contributes  only  the  zz  com¬ 
ponent  of  we  separate  it  by  writing 


s 

Azz 


(5.11) 


where  only  the  transverse  part  of  is  of  consequence  to  the  TE  problem  considered 
here. 

Using  (5.9)  in 


v.fi' =-Y'GV(£|£') 

we  easily  find  the  Lorentz  scalar  potential  Green’s  function  as 


ndK 


G<t(£lp')  =  ^  j 

n=l  n 


(5.12) 


(5.13) 


(which,  incidentally,  is  identical  with  G^„). 

The  eigenfunction  expansion  method  can  also  be  used  to  solve  (2.17)  for  the  Coulomb 
scalar  potential  Green’s  function,  with  the  result 

,.  dK 


n=l  ri 


(5.14) 
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The  integrations  over  the  radial  spectral  variable  k  in  (5.9)-(5.10)  and  in  (5.11)-(5.12) 
can  be  done  upon  using  [13] 


OO  j  1  /  \  *' 

f  Ju{kp)J^{kp')—  =  —  i  —  \  ,  i/>0 

J  K  /lU  V  j 

n  \f  / 


(5.15) 


r  // 

J  J^{kp)  Ju{kp')  =  -^J^{kp<)H]^\kp>),  i/>0 


and 


Trj 

2P 


MI=P<)H?\kp^) 


(5.16) 


iy>0  (5.17) 


where  /9>  denotes  the  greater,  and  /9<  the  lesser,  of  p  and  />',  and  is  the  Hankel  function 
•  of  the  second  kind  and  order  v.  In  the  above  and  henceforth,  we  omit  the  subscript  n  on  v 

for  notational  simplicity. 

Using  (5.17)  in  (5.10)  and  the  notation 


• 

c*  =  cos  y{'p  —  <p')  ±  cos  u{(p  +  ip') 

(5.18) 

—  sin  —  p')  ±  sin  u{p  +  p) 

(5.19) 

K'  =  7 

-^,JAkp<)H'-P(kp^)-va„ 

(5.20) 

• 

k 

bf  =  - 

■?pp' 

'  §-/Ukp')nl^Hkp]-^n„,  p>p' 

PMkp]H<^>'(kp'}  +  r(l„,  p<p' 
kp 

(5.21) 

• 

bf  =  < 

K  ' 

'  -^,Mf^p')Hl^^’{kp)-vnr,,  p>p' 

+  p<p 

(5.22) 

• 

(5.2.3) 

where  Jl  and  denote  derivatives  with  respect  to  the  arguments  of  the  respective  func¬ 
tions  and 
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fin  = 


P< 


Ttk'^pp'  \p>J 

we  can  express  the  elements  of  as 


(5.24) 


=  —  Y'  c~  y‘‘‘‘ 

(5.25) 

00 

(5.26) 

^  CO 

=  —  y  s-  b‘^^’ 

(5.27) 

oo 

(5.2S) 

Finally,  using  (5.15)  in  (5.14)  and  summing  the  resulting  series,  we  can  express  the 
Coulomb  scalar  potential  as 


Gi  =  -\nR+) 

27r 


where 


R^  =  (1  -2^  cos  (=^  +  4^)? 


(5.29) 


(5.30) 


with 


TT 


TT 


P<\i3 


0 


p>. 


(5.31) 


The  lamellar  part  of  can  also  be  transformed  by  a  procedure  analogous  to  that 
used  to  obtain  equations  (5.25)-(5.28).  When  the  corresponding  components  of  and 
are  added  to  form  the  transverse  part  of  the  Lorentz  vector  potential  Green’s  function, 
cancellations  occur,  resulting  in  the  relatively  simple  expressions 


TT 


G 


TT 


Ap<p' 


2;/? 


(5.32) 

(5.33) 
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(5.34) 


Gaw'  ~  t-  o  ^ 

VP  n=0 


where 


=  J  [j,Mi=p<)«i-i(kp>)  ±  J.Mkp<)«SAh>)] 

Finally,  using  (5.16)  in  (5.13),  we  express  the  Lorentz  scalar  potential  as 

_  OO 


(5.35) 


(5.36) 


(5.37) 


5.3  Series  Acceleration 

The  series  in  (5. 25)- (5. 28)  are  slowly  convergent,  which  can  be  appreciated  by  examining 
their  large-n  asymptotic  behavior.  Hence,  keeping  terms  up  to  0{n~^),  we  find 

b^/  ~  \t.  +  -f  2T+)1  (5.38) 


(5.39) 


bf  ^  (In  \T.  +  -iTl+2T+) 

bf  -  ±(ln  \t.  +  ^T.{T.  ±  2) 

6r'~±n„  [r_-{-^r_(r_T2) 
b'f  ~  -(In  [t_  +  ^(Ti  -  2T+) 

where  the  upper  and  lower  signs  correspond,  respectively,  to  p  >  p'  and  p  <  p',  and 

n.tel’.fe)’ 


(5.40) 


(5.41) 


Observe  that  the  leading  terms  in  (5.38)-(5.41)  are  0(n°).  The  series  (5.32)-(5.35)  and  (5.37) 
of  the  Lorentz  gauge  formulation  behave  slightly  better,  which  is  evident  from  the  asymptotic 
forms 


± 

(i„  ~  — 


i±  ^ 


Zi/TT  /  p< 


(5.43) 
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(5.44) 


Mkp<)  H'~^\kpy)  ~  tS- 

VK 

The  large-n  expansions  (5.38)-(5.41)  and  (5.43)-(5.44)  have  been  used  in  conjunction 
with  the  Kummer  transformation  [18]  to  accelerate  the  convergence  of  the  respective  series. 
In  this  method,  we  subtract  from  each  term  of  the  series  its 

large-n  form,  which  results  in  a  more  rapidly  converging  series.  To  preserve  the  correct 
value  of  the  sum,  the  modified  series  is  augmented  by  the  series  of  asymptotic  forms,  which 
can  be  summed  in  closed  form.  The  closed-form  sums  needed  in  our  case  are 

oo  cn 

— cos  nC^  =  — In (5.45) 

n=l  ^ 


E 


n=l 


r  . 

—  Sin  nQ 
n 


arctan 


(  sin  C* 
1  —  ^  cos 


(5.46) 


OO 

c  cos  nC^  = 

n=l 

^^sinnC*  = 

n=l 


^(cosC*  -0 

(  sin  C* 

~w~ 


(5.47) 

(5.48) 


These  sums  are  easily  derivable  from  the  geometric  series  [20]. 

After  the  Kummer  transformation,  all  series  have  terms  with  large-n  behavior  0(n~^), 
which  is  manageable.  However,  we  note  that  more  effort  is  required  to  accelerate  the  Coulomb 
gauge  series  than  is  needed  to  accelerate  the  Lorentz  gauge  series  to  the  same  rate  of  conver¬ 
gence.  On  the  other  hand,  the  scalar  potential  Green’s  function  of  the  Coulomb  formulation 
is  available  in  closed  form,  wherecis  that  of  the  Lorentz  formulation  is  in  a  series  form  and 
must  be  summed  numerically. 

Finally,  we  note  that  and  the  pp'  and  (pip'  terms  of  and  are  logarithmically 
singular  when  p  approaches  p'.  The  added  benefit  of  the  Kummer  transformation  is  the 
removal  of  the  singularities  from  the  corresponding  series. 


5.4  Incident  Field  and  Far  Field  Pattern 


The  TE  plane  wave  with  amplitude  Eq  impinging  upon  the  wedge  at  an  angle  ^pi  (Fig.  5.1) 
produces  in  tne  absence  of  the  arbitrary  scatterer  C  the  “incident”  field  The  com¬ 

ponents  of  E_\  which  appear  as  the  driving  function  in  the  integral  equation  (5.1),  are  easily 
found  as  [22] 

oo 

El  =  Eq-  e„e^(‘'+^^'/V'(A:y9)[cos  v{ip  -  pi)  +  cos  v{p  +  pi)]  (5.49) 

P  n=0 
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E'p  =  Eq-Y  ■r-e^^‘'+^''/V^(^^/9)[sin  v{p  -  pi)  +  sin  v{p  +  (/>.)] 

H  n=l 


(5.50) 


When  the  scatterer  is  introduced  near  the  wedge,  the  field  E',  is  modified  by  the  “scat¬ 
tered"  field  Ef ,  whose  source  is  the  surface  current  (of  density  J)  induced  on  C.  It  can 
be  shown  that  far  from  the  edge 

«  =  [^)  ’  (5-51) 

where 
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is  the  far  field  pattern,  in  which 
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Eip,  ^)  =  pY  rr-^n  -pY 

n=l 


n=0 


(5.52) 


(5.53) 


5.5  Sample  Results 

We  present  in  this  section  sample  numerical  results  for  several  simple  structures  near  a 
wedge  with  different  values  of  a  (see  Fig.  5.1).  In  all  cases,  /  =  300  MHz,  which  corresponds 
to  free-space  wavelength  A  =  1  m. 

In  Figs.  5.3(a)  and  5.3(b)  are  shown,  respectively,  the  current  distribution  and  the  far 
field  pattern  for  an  inverted  rectangular  trough  with  its  ends  touching  a  PC  ground  plane 
(which  corresponds  to  a  =  180°)  under  normal  incidence.  Our  results  are  compared  with 
that  obtained  for  the  trough  and  its  image  (a  rectangular  box)  by  a  free  space  code  [21], 


and  excellent  agreement  is  observed.  Since  the  Coulomb  and  Lorentz  gauge  based  programs 
produced  graphically  indistinguishable  results,  only  the  Lorentz  gauge  data  are  displayed  for 
this  and  other  cases  discussed  here. 

The  next  structure  considered  was  that  of  a  circular  cyhnder  of  radius  a  =  A/27r,  enclosing 
the  edge  of  a  wedge  with  a  =  30°,  thus  simulating  a  cylindrically  tipped  wedge,  for  which 
an  exact  solution  is  available  [24].  The  circular  profile  C  of  the  cylinder  was  approximated 
by  thirty  straight  line  segments.  The  current  distribution  on  C  for  (pi  =  45°  is  shown  in 
Fig.  5.4(a)  and  the  corresponding  far  field  pattern  in  Fig.  5.4(b),  where  the  exact  result  [26] 
is  also  shown  for  comparison.  When  the  wedge  angle  a  is  set  to  zero,  a  cylindrically  tipped 
half-plane  is  obtained.  The  corresponding  results  are  shown  in  Figs.  5.5(a)  and  5.5(b).  In 
this  case,  thirty  five  segments  were  used  to  model  the  cylindrical  profile  C. 

Finally,  in  Figs.  5.6(a)  and  5.6(b)  we  present  results  for  a  half-plane  (q  =  0°)  with  the 
edge  enclosed  by  a  square  box  with  the  diagonal  of  0.5A,  illuminated  by  a  plane  wave  with 
=  45°.  Since  the  edge  is  completely  shielded  from  the  incident  field,  its  position  inside  the 
box  should  have  no  effect  on  the  magnitude  of  the  induced  current.  This  is  indeed  observed 
in  Fig.  5.6(a),  where  one  result  is  for  a  half-plane  extending  to  the  center  of  the  box,  and 
the  other  is  for  the  case  where  the  edge  of  the  half-plane  is  displaced  O.lA  from  the  center 
of  the  box. 


n‘*  node 
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Figure  5.3:  (a)  Normalized  current  dens 
trough  on  a  ground  plane  illuminated  by 


8 


I  PM  I 


Figure  5.5:  (a)  Normalized  current  density  on  C  and  (b)  far  field  pattern  for  a  cylindrically 
tipped  half-plane  (a  =  0°)  illuminated  by  a  plane  wave  with  =  45°. 
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3.5 


Figure  5.6:  (a)  Magnitude  of  the  normalized  current  density  on  C  and  (b)  far  field  pattern 
for  a  half-plane  (q  =  0°)  with  the  edge  rncloserl  by  a  square  box  and  illuminated  by  a  plane 
wave  with  =  45°. 
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Chapter  6 
Conclusions 


In  this  report  the  Coulomb  gauge  formulation  and  solution  of  the  electromagnetic  field 
vector  and  scalar  potential  Green’s  functions  are  presented  for  both  waveguiding  and  scat¬ 
tering  structures.  For  the  scattering  cases  presented  here  the  Green’s  functions  are  placed  in 
mixed  potential  electric  field  integral  equations  (MPlE’s)  which  are  solved  by  the  moment 
method.  Green’s  functions  are  derived  for  free  space,  the  rectangular  waveguide  and  the  two- 
dimensional  perfectly  conducting  cylinder  and  wedge.  Comparisons  between  the  Coulomb 
and  Lorentz  gauge  formulations  are  made  in  terms  of  analytical  difficulty  and  computer 
computation  time  necessary  to  solve  for  the  current  on  scatterers  of  various  shapes. 

The  attractive  feature  of  the  Coulomb  gauge  is  that  for  the  cases  considered  here  the 
Coulomb  scalar  potential  is  either  mathematically  simplier  or  in  closed  form,  rather  than  in 
series  form  as  in  the  Lorentz  guage.  However,  this  advantage  of  the  Coulomb  gauge  is  to  some 
extent  offset  by  the  fact  that  the  resulting  vector  potential  series  converge  more  slowly  than 
the  corresponding  series  of  the  Lorentz  gauge.  Therefore  acceleration  of  these  series  requires 
more  analytical  work  than  is  necessary  with  the  corresponding  Lorentz  gauge.  Also,  even 
in  cases  where  the  Lorentz  vector  potential  Green’s  dyadic  is  diagonal,  its  Coulomb  gauge 
counterpart  is  likely  to  be  fully  populated.  This  is  true,  for  instance,  for  both  free  space 
and  the  rectangular  waveguide.  However,  for  the  cylinder  and  the  wedge  both  gauges  have 
fully  populated  vector  potential  Green’s  functions.  Consequently,  for  these  two  ca^es,  the 
Coulomb  gauge,  MPIE’s  are  noticeably  more  efficient  than  their  Lorentz  gauge  counterpart, 
since  they  do  not  require  summation  of  scalar  potential  series.  For  example,  for  both  the 
cylinder  and  the  wedge,  when  -fo  basis  functions  are  used  to  represent  the  current,  a  15% 


savings  in  matrix  fill  time  was  realized  when  the  Coulomb  gauge  was  used  instead  of  the 
Lorentz  gauge. 

To  summarize,  this  research  has  shown  that  the  Coulomb  gauge  may  be  perferable  to  the 
Lorentz  guage  when  the  MPIE  formulation  is  used  to  solve  scattering  problems  in  complex 
environments  which  admit  closed  form  solutions  of  the  Poisson  equation.  The  number  of 
such  geometries  is  admittedly  limited.  In  addition  to  the  cases  presented  here  the  only  other 
three  dimensional  problem  known  to  the  authors  which  falls  into  this  catagory  is  that  of  the 
scatterer  or  antenna  of  arbitrary  shape  in  the  presence  of  a  sphere. 

It  should  be  pointed  out,  however,  that  the  Coulomb  scalar  potential  contains  the  static 
limit  portion  of  the  electromagnetic  field.  One  can  therefore  isolate  and  observe  the  physical 
behavior  of  that  part  of  the  field  (the  Coulomb  vector  potential  portion)  that  approaches 
zero  as  the  excitation  frequency  is  decreased.  This  observation  may  prove  useful  in  electro¬ 
magnetic  compatability  studies. 

In  conclusion,  a  major  benefit  of  this  research  is  that  we  have  been  able  to  demon¬ 
strate  tliat  it  is  not  necessary  to  solve  classical  electromagnetic  scattering  problems  with 
the  Lorentz  gauge.  The  Coulomb  gauge  has  been  shown  to  be  one  alternative,  although 
not  genei  ally  an  improvement.  It  has  recently  been  suggested  [25]  that  the  Coulomb  gauge 
might  reduce  numerical  instabilities  in  time  domain  integral  equation  scattering  analysis. 
It  has  also  been  noted  by  the  authors  that  in  the  physics  literature  (e.g.  [2])  the  Coulomb 
gauge  is  applied  in  momentum  space  which  is  analogous  to  the  spectral  domain  in  classical 
electromagnetism,  d  ime  domain  and  spectral  domain  applications  of  the  Coulomb  gauge 
remain  possible  avenues  for  future  research. 
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